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Abstract —A new class of folded subspace codes for non¬ 
coherent network coding is presented. The codes can correct 
insertions and deletions beyond the unique decoding radius 
for any code rate R € [0,1]. An efficient interpolation-based 
decoding algorithm for this code construction is given which 
allows to correct Insertions and deletions up to the normalized 
radius s (1 — ((l//i + h)/{h — s + 1))R), where h is the folding 
parameter and s < h is a decoding parameter. The algorithm 
serves as a list decoder or as a probabilistic unique decoder 
that outputs a unique solution with high probability. An upper 
bound on the average list size of (folded) subspace codes and on 
the decoding failure probability is derived. A major benefit of the 
decoding scheme is that it enables probabilistic unique decoding 
up to the list decoding radius. 

Index Terms —Network coding, subspace codes, lifted MRD 
codes, folded subspace codes 

I. Introduction 

Subspace codes have been proposed for error control for 
noncoherent random linear network coding, e.g. when the 
network topology and the in-network linear combinations are 
not known by the transmitter and the receiver [Ij, [2j. Kotter 
and Kschischang proposed a Reed-Solomon like construction 
based on rank-metric codes (referred to as KK codes) that can 
be decoded efficiently [IJ. List decodable variants of subspace 
codes have been proposed in [3j-[9] and allow to correct 
insertions and deletions beyond half the minimum subspace 
distance. The challenge of list decoding subspace codes is to 
decrease the size of the list of candidate codewords, which is 
exponential in the dimension of the transmitted subspace [ 1 ^. 
Most list decodable subspace codes are based on KK codes and 
control the size of the list by restricting the message symbols 
or the code locators to belong to a subfield. Guruswami and 
Wang [9] showed that punctured subspace codes can be list 
decoded up to the theoretical limit for any code rate. The list 
size for this decoder is further reduced in by applying 
hierarchical subspace evasive sets. The output of this decoder 
is a basis for the affine space of candidate solutions resulting 
in a very large list of exponential size in the dimension of the 
transmitted subspace with high probability. 

In this paper we define a new class of folded subspace codes 
that can be decoded from insertions and deletions for any code 
rate. We present an interpolation-based decoding algorithm 
that can be used as a list decoder or as a probabilistic unique 
decoder. Both schemes can correct insertions and deletions 
beyond half the minimum subspace distance for any code rate. 


The probabilistic unique decoder returns a unique solution 
with high probability and requires at most 0{s^nf) operations 
in ¥qm, where s is a decoding parameter (small integer) and 
rir is the dimension of the received subspace. The decoding 
scheme is well suited for practical application. We give an 
upper bound on the probability of a decoding failure (i.e. a list 
of size larger than one) and verify the results by simulations. 

This paper is structured as follows. In Section II, we 
describe the notation and give basic definitions. Section III 
introduces a new class of folded subspace codes and presents 
an efficient interpolation-based decoding scheme. In Sec¬ 
tion IV we apply the algorithm to list and unique decoding 
of folded subspace codes and highlight the improvements of 
the decoding scheme. Finally, Section V concludes this paper. 

II. Preliminaries 
A. Finite Fields and Subspaces 

Let q be a power of a prime, and let F, be the finite field 
of order q and let F^m be its extension field of degree m. 
Any element from F,™ can be represented by a row vector 
of length m over F, for a fixed basis. By F^ we denote 
a vector space of dimension N over F^ and the set of all 
subspaces of F^ is the projective space 'Pq{N). The set of all 
^-dimensional subspaces of F^ is the Grassmannian and is 
denoted by 0q{N,i). We denote matrices and vectors by bold 
uppercase and lowercase letters such as A and a and index 
their elements beginning from zero. The rank of a matrix A € 
■pmxn denoted by rk(A) and the kernel of A is denoted 
by ker(A). The row space of a set of vectors B over F, is 
denoted by (B)^. For two subspaces U,V G Vq{N), the direct 
sum U ®V is the smallest subspace containing both U and V. 
The subspace distance between IA,V in Vq{N) is 

dsiJA, V) = dim(^^) -f dim(V) — 2 AimiU fl V). (1) 

A subspace code is a nonempty subset of Vq{N), and has 
minimum subspace distance dg when all subspaces in the code 
have distance larger than or equal to dg from each other. 

As channel model we use the operator channel from [IJ. 
Such a channel has input and output alphabet Vq{N). The 
output U is related to the input V with dim(V) = rit by 

U = nn,-5{V)®£ ( 2 ) 

where FLnt-siV) returns a random (nj — (5)-dimensional sub¬ 
space of V, and £ denotes an error space of dimension 7 with 


V n f = 0. The distribution of 'Hn^-s{V) does not affect the 
performance of the code and can be chosen to be uniform 
(see [IJ). The dimension of the received subspace U is thus 
Ur = nt —6 +j and we call 6 the number of deletions and 7 
the number of insertions. 

B. Linearized Polynomials 

For any element a € F^m and any integer i let al'l a‘^' 
be the Frobenius power of a. A nonzero polynomial of the 
form p{x) = Pi ^ Pd ^ 0 , is called a 

linearized polynomial of 5 -degree deg^(p(x)) = d, see illl . 
[12]. Evaluating a linearized polynomial forms a linear map 
over Fq, i.e. for all a, 6 e F^ and a:i,a ;2 € F^m, we have 
p{axi -f bx 2 ) = ap{xi) + bp{x 2 ). The noncommutative com¬ 
position f(x)®g{x) = f{g{x)) of two linearized polynomials 
f{x) and g(x) of 5 -degree di and d 2 is a linearized polynomial 
of 5 -degree di+d 2 . The set of all linearized polynomials over 
¥gm forms a noncommutative ring Lqm[a;] with identity under 
addition “H-” and composition The Moore matrix of the 


vector a = (oq oi 

• ■ O^n — l 

G F”m is defined as 



( ao 

ai 

^n —1 \ 


M,.(a) = 

„W 

aW ._ 


(3) 


W-‘' 

[r'-l] 

a[ 

■ «t;'/ 


The rank of (a) 

is min{r 

n} if the elements oq 

,..., dji—i 


are linearly independent over F^, see [ 121 . 

III. Interpolation-Based Decoding oe Folded 
Subspace Codes 

We present a new construction of folded subspace (FS) 
codes that can be decoded from insertions and deletions 
beyond the unique decoding radius for any code rate R. This 
work is motivated by the constructions in [131 and [141 . 

Let a be a primitive element of the field F^m with polyno¬ 
mial basis a°, a^,..., over Fg. 


Lemma 1 The minimum subspace distance of the code 
FSub[/i; rit, fc] is ds.min = 2 (nt-r!l+l)- 


Proof: Let V and V' be two distinct codewords generated 
by f{x) and g{x) with 5 -degrees less than k and suppose 
dim(V n V') > [-^l. Then f{x) and g{x) must agree on 
> k linearly independent points, which is not possible 
since the 5 -degree of both polynomials is less than k. Thus 
the dimension of the intersection space V Cl V' can be at most 
[|] — 1. Using (1) we have 


d,{V, V') = 2nt-2 dim(V n V') = 2{nt - 


1 ). 


A. Interpolation Step 

Suppose we receive a basis of dimension Ur = rit — 6 + j 
{{xj,yj,o,Vj.i,■■■,Vj.h-i) ■■ j e [ 0 ,nr - 1 ]} 
of the received subspace U. Let the matrix 
[x^, ..., G contain this basis 

as rows. Suppose we receive a F^-linear combination of the 
transmitted basis vectors of the form 

rtt —1 

^ A, (a^■^ ,..., 

j=o 

with Xj G Fq. Due to the linear property of linearized 
polynomials we can rewrite this as 

C Tit — l rtt — l rtt — l \ 

J2x,a^\f{ ^A,. (4) 

i=o j=o j=o ) 

For the interpolation step we must solve the following 
problem. 

Problem 1 Given the integers D and s < h,l < s < h, find 
a nonzero (s + 1)-variate linearized polynomial of the form 


Definition 1 An h-folded subspace code FSub[h; nt, fc] of 
dimension rit, where hrit < m, is defined as the set of 
subspaces 

..., : j G j})^ 

for all f{x) G ¥,qm[x\, degg{f{x)) < k, where J = [0,nt —1]. 
The dimension of the ambient space 

1U,= 0.-.0Fgm 

h times 

is iV = nt 0 hm, since the vectors in the space 
..., have nonzero components at the nt 

known positions 0, h, 2h ,..., {rit — l)h only. The zeroes at 
the known positions do not need to be transmitted and can be 
inserted at the receiver. The code rate is R= — , —r. 


Q{x,yi,...,ys)= Qq{x) 0 Qi{yi) H-h Qs{ys), (5) 

which satisfies for all i G [0, /i — s], j G [0, rir — 1].' 

• (f{xj(y.\yj^i^yjyj^\^... ^yj^ij^s—i) = 0 , 

• degq(Qo(a;)) < D, 

• degq(Q(!(yf)) < £» - (/c - 1), G [1, s]. 

Here we use (0) to determine the code locators for the {h — 
s 0 1) interpolation tuples for each dimension as XjO^^'di G 
[0,h — s],j G [0,711- — 1]. A solution to Problem X can be 
found by solving a homogeneous linear system of equations. 
Denote the polynomials of (5) by Qoix) = and 

QiiVi) = Let the matrix T contain all rirfh — 

s 0 1) interpolation tuples {xja\ yj^i, yj,i+i ,..., yj^i+s-i). 
Mi G [0, h — s], j G [0, Ur — 1] as rows and denote by the 
£-th column of T for i G [0,s]. The coefficients qij can be 
found by solving a linear system 

R • qf = 0 


(6) 








where R is an nr{h — s + 1) x D{s + 1) — s{k — 1) matrix; 

R= (M,3(t^)^,Mi3_fe+i(tf )^,..., M,5-fc+i(tJ)^) (7) 
and q/= (qo,o, ■ ■ ■, go,D-i| \qs,o, ■ ■ ■ ,qs,D-k)- 


Lemma 2 A nonzero polynomial fulfilling the interpolation 
constraints in Problem 1_ exists if 

nr.(/i — s + 1) + s(/c — 1) + 1 


Proof: Problem X forms a homogeneous linear system 
of nr{h — s + l) equations in Z?(s + 1 ) — s(fc— 1 ) unknowns. 
This system has a nonzero solution if the number of linear 
independent equations is less than the number of unknowns, 
i.e., if 

nr{h — s + 1) < D{s + 1) — s{k — 1) (9) 

^ ^ nr{h — s + 1 ) + s{k - 1 ) + 1 

s + 1 ^ 

The receiver knows the code parameter k and the 
decoding parameter s and can compute the degree restriction 
D in ( 8 ). 


Theorem 1 Let Q (x,yi,..., ys) ^ 0 fulfill the interpolation 
constraints in Problem L If 

then 

P(x) =*' Q{x,f{x),f{ax),...,f{a'‘~^x))=0. (11) 

Proof: The dimension of the noncorrupted subspace Vr\U 
is rit — S. The code locators of the noncorrupted dimensions 
are linearly independent and thus we have (rit — 5){h — s + 1) 
linearly independent interpolation points (roots) in (^. Since 
deg^(P(a;)) < D the dimension of the root space of P{x) is 
at most 79 — 1. If 

D < [rit — 5)[h — s + 1) (12) 

then P[x) has more linearly independent roots than its degree. 
This is possible only if P[x) = 0. Combining (9) and (12) 
and using 

Ur =nt+^ — 5 (13) 


we get ( 10 ) . ■ 

For h = s = 1 (no folding) we have y+6 < rit — k+l which 
is identical to the decoding radius of KK codes in [IJ. For s = 
h the algorithm in [131 is identical to the proposed scheme. 
By using the approximation R ss nt(nt+^) normalized 
decoding radius Tf = is given by 


V 


s 


nt + hm \ 
m[h — s + 1) ) 


If nth « m then we may write t/ Ri s ^1 — jpz^^R^ ■ 
Problem X can be solved by the efficient interpolation 
algorithm in [151 requiring at most 0{s^nrD[h — s + 1)) < 
0[s'^n^) operations in F^m. 


B. Root-Finding Step 

Given a polynomial Q [x,yi,... ,ys), we must find all 
polynomials f[x) G Lqm[x] of degree less than k which are a 
solution to (11) . To increase the probability to find a unique 
solution we use a similar idea as in [Q, [ 8 J. The solution space 
of the interpolation system (© has dimension larger than one in 
general. In this case, there exists a set of linearly independent 
linearized polynomials Q {x,yi,... ,ys) which are a solution 
to Problem X Instead of one polynomial we use a basis for 
the solution space of ( 6 ) to increase the probability that the 
root-finding system has a unique solution. We now derive a 
lower bound on the dimension of the solution space of ( 6 ). 

Lemma 3 The dimension d/ of the solution space of the inter¬ 
polation system ( 6 ) satisfies dj > s[D — /c + l) — 7 (/i — s + 1). 

Proof: Let T' contain the [nt — 5){h — s -\- 1) non¬ 
corrupted interpolation tuples as rows and denote by the 
7-th column of T' for i G [0,s]. Assume w.l.o.g. that the 
first [nt — 6)[h — s -G 1) rows of R correspond to the 
noncorrupted interpolation tuples and denote this matrix by 
R'. The first D columns of R' form a [nt — S)[h — s-G1) x D 
Moore matrix M£)(tQ^)^ of rank D since the elements in 
tg are linearly independent and (12) holds. The s Moore 
matrices M£,_(fc_i)(t'i ^)^,... ,M.D_(k-i)[^'s'^Y' are linear 
combinations of the rows of M£)(tQ^)^ and hence do not 
increase the rank. (/i — s + 1 ) interpolation constraints (i.e. 
rows) are added to R' for every malicious dimension. Thus 7 
insertions can increase the rank of R' by at most y[h — s + 1). 
Hence we have rk(R) < D-Gy[h—s-Gf). The dimension of the 
solution space of the interpolation system di :=dimker(R) is 

di > 79(s+l) —s(fc —1) —rk(R) = s(79 —fc+1)— 7 (/i—s + 1). 


We now set up the root-finding system using dj polynomials 
Q [x,yi,... ,ys). Define the polynomials 


sf^x) = q^t + q^^^Jx + qf^x"^ + • ■ • + qfjx^^ 


for 


[l,d/] and the vectors 


[a^^)... and qo,i = ^ 

for i, j G [0, k — 1], The root-finding matrix is 


%,i 


/ b, 


B = 


0,0 


Ki b 


,[-i] 

1,1 

J-(fe-i)] 


[- 1 ] 


0,1 

-(fe-i)] 


\"fc-l,fe-l '^k-2,k-l ■ • ■ / 




(14) 


anHn-Z'n 

and q — (^qo,o Bo,! ■ • ■ Bo.fe-i J ■ 

We can find the coefficients of the message polynomial f[x) 
by solving the linear system 


Bf = -q (15) 

where f[x) is connected with the vector f by f = 

^/o f[ ... flJ'i . The root-finding system (15) has 
at least one solution, i.e. q is always in the column space of B 













since we guarantee that the transmitted message polynomial 
f{x) is a solution to (111 if 7 and 6 satisfy (101 . Due to the 
lower triangular structure of B the root-hnding system (151 
can be solved in at most 0{k‘^) operations in F^m. 

IV. List and Unique Decoding oe Folded Subspace 
Codes 

We now show how the interpolation-based decoding scheme 
from Section III can be used as a list decoder and as a 
probabilistic unique decoder. We focus on the probabilistic 
unique decoding approach that is well suited for applications. 

A. List Decoding Approach 

The solution space of the root-hnding system (151 is an 
affine subspace over Fq. In case B in (151 has rank less than 
k, we obtain a list of possible message polynomials f{x) that 
satisfy ( 111 . 

Lemma 4 The dimension of the affine solution space of (151 
is at most 


Ur + rit in (171 we obtain the limits ui = and Um = 

j ^ 

Theorem 2 Let FSub[/i; nt, fc] be a constant dimension sub¬ 
space code over F^m and let N = rit hm be the dimension 
of the ambient vector space. Let the number of insertions 7 
and deletions 5 fulfill (101 . The average list size Lf^r), i.e. 
the average number of codewords at subspace distance at most 
T = 7 -f S(5 from a received Ur-dimensional subspace satisfies 

L{t) < 1 -f 

Proof: Let the received subspace y be chosen uniformly 
at random from all subspaces in the Grassmannian GqiN, Ur). 
The number of n*-dimensional subspaces in subspace distance 
at most r from y is Vs{nr, rit, t). If r satishes (101 we know 
that the causal (transmitted) codeword is in subspace distance 
at most T from y. There are — 1 noncausal codewords 
(subspaces) out of [^] possible nt-dimensional subspaces. 

Thus there are on average 


Proof: The lower triangular root-hnding matrix B has full 
rank if and only if all diagonal elements bop,..., bo,fe_i are 
nonzero vectors. The entries of each bp^ are the evaluations of 
di polynomials of degree at most s — 1 at 1 € [0, fc — 1]. 
Since the conjugates a,a^^\ ... are all distinct and 

deg(i7o^^(a:)) < s for all £ G we can have bo,i = 

0,i € [0, fc — 1] at most (s — 1) times. For each bo,i = 0,i € 
[0,k — 1] the coefficient can be any element in F^m. Thus 
the dimension of the affine solution space is at most 

■ 

Using di polynomials for the root-hnding step does not 
reduce the worst case list size. The probability that B is 
nonsingular increases with dj and thus the average list size is 
reduced. 


Lemma 5 The number Vs{nr,nt,T) of nt-dimensional sub¬ 
spaces in 'Pq{N) at subspace distance at most t from a fixed 
Ur-dimensional subspace in Vq{N) is 


Vs{nr,nt,T) = ^ 


Tlr-Vrit) 

Tlr 

N — nr 


J . 

j - Hr + rit_ 


(16) 


where ui = and Um = 


Proof: Denote by Ns{nr,nt,t) the number of nt- 
dimensional subspaces in subspace distance exactly t from an 
nj.-dimensional subspace. In [1^ Lemma 2] it is shown that 

iV 5 (n„nt,f)=g“(*-“)H[^“H (17) 


if n = is an integer and 0 otherwise. The number of 

all nt-dimensional subspaces at distance at most r is 


Vs{nr:nt,T) = y^Ns{nr,nt,t). 

Since Ns(nr,nt,t) 7 ^ 0 if and only if n = is an 

integer we rewrite (171 in terms of u. Substituting t = 2u — 


L'{t) = 

LntJ 

noncausal codewords in subspace distance at most r from the 
received subspace. Let ui = ] and Um = 

Using Lemma ^ and the approximation < [;] < 

4g^(” (see [IJ) we have 


l'{t)= {^qmk _^:ys{nr,nt,T) 


r] 

LntJ 


r,mk 


< 


qnt(N—nt) 




.j{j-nr+nt) 


Ur 

N — nr 

.3 _ 

J -rir-\- nt_ 




<i + l 


<q 


,mk—nt{N-nt) 


(U‘) 


N — Tlr 

Um 

{Um — nr + nt) 


Tlr 


Um_ 



=16 
= 16 
= 16 


, /^qUminr—Um) . —nr-l-nt)(N —Hr — (tlm —"r-l-nt)) 

2^^ qmk—nt(N—nt)-\-Umnt-\-{,Um—nr+nt)(N—Um—nt) 

qmk+{nr-Um){nt+Um-N) 




mfc-r(rar—[ "*' J)("* + L~ 


:j-A) 


Including the causal codeword we get L{t) = 1-\-L (r). 


B. Probabilistic Unique Decoder 

In the worst case the decoder outputs an exponential number 
of candidate message polynomials. We show that a list of size 
larger than one is a rare event. This allows us to use the 
algorithm as a probabilistic unique decoder which returns a 
unique solution or a decoding failure in case the list size is 
larger than one, i.e. rk(B) < k. 


































The root-finding system (15) has a unique solution if the 
rank of B is full. This is fulfilled if and only if at least one 
entry of each bg.i, i & [0, fc — 1] is nonzero. 


Lemma 6 Denote by dj the dimension of the solution space 
of (^. Then the decoding failure probability is upper bounded 
by 


Pe < k 





( 18 ) 


under the assumption that the coefficients of the polyno¬ 
mials Bq^\x),..., are independent and uniformly 

distributed over F^m. 


Proof: Evaluating Bq^\x),£ € at the distinct 

elements a, , • ■ •, gives a codeword of a (fc, s) Reed- 

Solomon code Crs- The probability to get a unique solution 
is then equal to the probability to get a weight k codeword. 
Since similar to and [8] we assume that the coefficients of 
B^\x),£ G are independent and uniformly distributed 

over Fqm. Hence we get a uniform distribution over the code 
book of Cfis- Using the approximation from [17. Equation 1] 
the probability Ps to get a codeword of full weight k is 

1 


no. of vectors of weight k in F^,, 
total no. of vectors in F^m 


= 1 



The probability that one in bo,i is zero is at most 

1 — Ps. The probability that one bo,i = 0,z G [0,fc — 1] is 
upper bounded by 


Pr[bo,. = 0] < (l-Ps)-^^ = 




The probability that at least one bo,i = 0 for i = 0,..., fc — 1 
is thus upper bounded by 


rfc-i 


Pe < Pr 


U bo,. = 0 

.i=0 



We restrict dj to be larger than a threshold p, i.e. p < dj, 
and get 


Ds > 'yfh — s -I- 1) -f s{k — 1) p. (19) 


To ensure that f{x) is a root of P{x) in (11) the degree D 
must satisfy (12) . i.e. D < {rit — 6){h — s 1). By combining 
(12) and (19) we get 


s{nt{h - s1) - {k - 1)) - p 

j + sd < ---—- 

h — s 1 


( 20 ) 


Under the assumption that the coefficients of the polynomials 
B^\x),£ G are independent and uniformly distributed 

over Fqm we can use Lemma ,6 to upper bound the failure 
probability. If 7 and 5 fulfill (20) we can find a unique solution 
f{x) satisfying (11) with probability at least 



( k 1 ^ TTh 

Using R Ri ntint-i-hm) normalized decoding radius r„ = 
of the probabilistic unique decoding approach is 

nt + hm \ _ p 

m{h — s1) ) {h — s-\-l)nt 

In a setup where nth Ri m we have t„ < s(l — ~ 

p/{,{h - s -I- l)nt). 

To adjust the decoding radius at the receiver we express the 
degree constraint D in terms of p. Combining (12) and (19) 
we get nr(/i —s-l-l)-l-s(A: —l)-|-/i < (s-|-l)(nt ——s-l- 1). 
Erom (12) we get D{s + 1) < (s -f l){nt — 6){h — s 1) and 
choose 


< s 1 - 


D = 


Urfh — s -f 1) -f sfk — 1) p 

s + l 


The computational complexity of the unique decoder is domi¬ 
nated by the interpolation step, which can be solved requiring 
at most 0{s^nrD{h — s -f 1)) < O(s^n^) operations in F^m 
using the efficient algorithm in [^. The pseudo code for the 
probabilistic unique decoder is given in Algorithm X 


Algorithm 1: UniqueDecodeFS(x^, ..., 

Input : A basis (x^, ..., for the 

n^-dimensional received subspace 
Output: A polynomial f{x) G Lgm[a;] : deg(/(a;)) < k 
or “decoding failure” 

1 Set up T G to contain all interpolation 

tuples of Problem X as rows and denote by to,..., tg the 
columns of T 

2 Interpolation step: 

3 ... ,(5^“^^^ InterpolateBasis(t^, tf’,..., t^) 

Root-finding step: 

4 Q* = {QW : deg,(gW) < G [l,d/]} 

5 Set up the root-finding matrix B as in (14) using all 
polynomials in Q* 

6 if bo.i 7 ^ 0, Vi G [0, k — 1 ] then 

7 Solve B • f = qo and define f{x) from f 

8 Output: f{x) 

9 else 

10 1^ Output: “decoding failure” 


C. Performance Analysis 

We compare the performance of our proposed code con¬ 
struction with the code constructions by Kotter and Kschis- 
chang [IJ, Mahdavifar-Vardy [131 and Guruswami-Xing [9J. 
For a fair comparison we select the code parameters such that 
each codeword contains the same number of symbols. Figure X 
shows that the code by Mahdavifar and Vardy only can correct 
errors for very small rates. The construction by Guruswami 
and Xing achieves the best decoding radius for all rates but 
puts out a very large list with high probability. 


1 - k 


















The proposed code construction can correct insertions and 
deletions for all code rates and returns a unique solution 
with high probability, which is a major benefit for practical 
applications. 
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Fig. 1. 



D. Simulation Results 

Consider a folded subspace code with parameters rit = 3, 
h = 3, k = A, s = 2, q = 2 and m = hrit = 9. For p = 1 (maximum 
decoding radius) we simulated 5.2 • 10® transmissions over an 
operator channel with <5 = 0 deletions and 7 = 2 insertions 
and observed a fraction of 7.80 • 10“^ decoding errors (upper 
bound 3.1310“^). For /r = 2 the code corrects 5 = 0 deletions 
and 7 = 2 insertions. We simulated 5.5-10® transmissions over 
an operator channel with parameters 5 = 0 , 7 = 2 for /r = 2 and 
observed a fraction of 1.97-10“® decoding errors (upper bound 
2.44-10“"^). For /i = 3 the code can correct 5 = 0 deletions and 
7=1 insertions. After simulating 1.62 • 10^ transmissions for 
7 = 1 we observed no decoding failure so far (upper bound 
1.91-10“®). The simulation is still in progress. The simulation 
results show that the assumptions in Lemma 6_ are reasonable. 

The decoder in [13] can not correct any insertions and 
deletions for these code parameters, if the same number of 
symbols is transmitted. The code from can correct the same 
number of insertions and deletions for the given parameters 
but will output a large list with high probability instead of a 
unique solution. 

V. Conclusion 

A new family of folded subspace codes for error correc¬ 
tion in noncoherent network coding scenarios was presented. 
The codes are more resilient against injections of malicious 
packets (insertions). An efficient interpolation-based decoding 
algorithm was presented that can be used as a list decoder 
or as a probabilistic unique decoder. The decoder corrects 
insertions and deletions for any code rate. We showed that the 
probabilistic unique decoder outputs a unique solution with 
high probability and gave an upper bound on the average list 
size and on the decoding failure probability. The decoding 


radius of the unique decoder can be adjusted to control the 
decoding radius vs. failure probability tradeoff. 
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